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ABSTRACT 



Context. We present the analysis of Ho-3, an Ha narrow-band imaging follow-up survey of 409 galaxies selected from the HI Arecibo 
Legacy Fast ALFA Survey (ALFALFA) in the Local Supercluster, including the Virgo cluster, in the region 1 l h < RA. < 16 h ; 4° < 
Dec. < 16°; 350 < cz < 2000 km s -1 . 

Aims. Taking advantage of Ho-3, which provides the complete census of the recent massive star formation rate (SFR) in Hi-rich 
galaxies in the local Universe and of ancillary optical data from SDSS we explore the relations between the stellar mass, the HI 
mass, and the current, massive SFR of nearby galaxies in the Virgo cluster. We compare these with those of isolated galaxies in the 
Local Supercluster, and we investigate the role of the environment in shaping the star formation properties of galaxies at the present 
. cosmological epoch. 

Methods. By using the Ho- hydrogen recombination line as a tracer of recent star formation, we investigated the relationships between 
atomic neutral gas and newly formed stars in different environments (cluster and field), for many morphological types (spirals and 
dwarfs), and over a wide range of stellar masses ( 10 7 5 to 10' 1 5 M ). To quantify the degree of environmental perturbation, we adopted 
an updated calibration of the HI deficiency parameter which we used to divide the sample into three classes: unperturbed galaxies 
^ ■ (Def H i < 0.3), perturbed galaxies (0.3 < Def HI < 0.9), and highly perturbed galaxies (Def HI > 0.9). 

Results. Studying the mean properties of late-type galaxies in the Local Supercluster, we find that galaxies in increasing dense local 
galaxy conditions (or decreasing projected angular separation from M87) show a significant decrease in the HI content and in the 
mean specific SFR, along with a progressive reddening of their stellar populations. The gradual quenching of the star formation 
occurs outside-in, consistently with the predictions of the ram pressure model. Once considered as a whole, the Virgo cluster is 
' effective in removing neutral hydrogen from galaxies, and this perturbation is strong enough to appreciably reduce the SFR of its 

( entire galaxy population. 

Conclusions. An estimate of the present infall rate of 300-400 galaxies per Gyr in the Virgo cluster is obtained from the number of 
existing Hi-rich late-type systems, assuming 200-300 Myr as the time scale for HI ablation. If the infall process has been acting at a 
constant rate, this would imply that the Virgo cluster has formed approximately 2 Gyr ago, consistently with the idea that Virgo is in 
a young state of dynamical evolution. 
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1. Introduction nearby galaxies. This observational effort is complemented by 

new progress in numerical simulations and theory in modeling 

A deeper understanding of the basic processes that regulate the the proC esses that regulate the formation of new stars in galaxies 

life of galaxies, i.e, the conversion of gas into new stars and ( e g Somerville & Primack 1999) 

their progressive aging, is becoming a reality owing to the recent ALFALFA (Giovanelli et al. 2005) is a blind survey that just 

multi-wavelength surveys that are providing a wealth of infor- ended at Arecibo (October 2012), aimed at obtaining a census 

mation on the different phases of the interstellar medium (ISM), of HI sources within 70 00 sq degrees of the sky accessible from 

the stellar content, and the current star formation rate (SFR) in Arecibo (including the Virgo cluster), with a sensitivity of about 

^3 mJy km s _1 , corresponding to 10 7 7 M of HI at the distance 

* Based on observations taken at the observatory of San Pedro Martir of Virgo. A catalog listing 40% of the whole ALFALFA sources 

(Baja California, Mexico), belonging to the Mexican Observatorio 

Astronomico Nacional. 1 2 mly is the noise per channel after Hanning smoothing to 10 

** Hubble Fellow km s~' (see Haynes et al. 201 1). 
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Fig. 1. Sky distribution of 409 Hl-selected galaxies in Ha-3, coded according to their morphology (E-SO-SOa: circles; giant spirals 
Sa-Sm: squares; Irr-BCD: triangles) The two vertical dashed lines mark the adopted boundaries of the Virgo cluster. Red symbols 
refer to galaxies considered as part of the Virgo cluster, blue are galaxies considered as isolated. Superposed are the X-ray contours 
from ROS AT (Bohringer etal. 1994). 



(a.40 catalog) has been published by Haynes et al. (2011). The 
Sloan Digital Sky Survey (SDSS, York et al. 2000) revolution- 
ized our knowledge on the stellar content of galaxies in the lo- 
cal Universe. The GALEX mission (Martin et al. 2005) with its 
all-sky survey (and other smaller but deeper surveys, such as 
GUVICS by Boselli et al. 2011) disclosed our view of the ul- 
traviolet sky, thereby providing quantitative knowledge on the 
massive, recent star formation in galaxies. 
Based on these datasets Huang et al. (2012) were able to deter- 
mine the scaling relations between the atomic hydrogen content, 
the stellar content, and the SFR in HI selected galaxies over a 
range of mass approaching one million (from 10 6 to 1O 12 M ). 
Similar analyses from optically selected surveys were carryed 
out in the local volume (Lee et al. 2007, 2009, Bothwell et 
al. 2009) and beyond (e.g. GASS, the GALEX Arecibo SDSS 
Survey, Catinella et al. 2010; Schiminovich et al. 2010), limited 
to galaxies with stellar masses in excess of 10 10 M however. 
Huang et al. (2012) based their study on the star formation 
rate estimated by fitting spectral energy distributions (SEDs) ob- 
tained by combining optical (SDSS) and UV (GALEX) photom- 
etry. This technique requires that extinction in the UV is prop- 
erly accounted for (Cortese et al. 2006, 2008b) and depends on 
the assumed star formation history (SFH) and metallicity. On 
the other hand the UV or Ha luminosity alone (see Kennicutt 
1998) provides a more istantaneous estimate of the SFR and re- 
quires that the SFH of galaxies has remained stationary over the 
time scale typical of A (10 8 yr) or OB stars (10 7 yr). This ap- 
proach is more suitable when it is necessary to determine the 
SFR of galaxies that are subject to abrupt changes in their SFH, 
e.g., when they suffer from environmental disturbances, such as 
quenching by ram pressure. Vollmer et al. (2004) and Boselli et 
al. (2006) have shown that in massive galaxies, such as NGC 
4569, ram pressure can produce complete quenching of the star 
formation on a timespan as short as 100-300 Myrs, and even 
shorter in dwarf galaxies (Boselli et al. 2008). 
The SFR determined from UV magnitudes has been compared 
to that from recombination lines by Boselli et al. (2009), who 
arrived at the conclusion that the accuracy of SFR based on 
UV (GALEX) luminosity critically depends on the quality of 
the data used, i.e, whether they come from the shallow (100 sec 



exposure) All-sky-survey (AIS) or from the deeper MIS survey 
(typically 1500 sec integration). Huang et al. (2012) reported 
that 25% of the a A targets do not have a counterpart in GALEX, 
52.2% of the remaining match to UV sources found in the AIS, 
and only 23.1% matches those in the MIS. This severely ham- 
pers an unbiased comparison of the star formation properties of 
galaxies in the Virgo cluster with respect to isolated galaxies 
in the Local Supercluster, i.e, the main goal of the present in- 
vestigation using UV data, because the large majority of Virgo 
galaxies have been observed in the MIS (owing to GUViCS by 
Boselli et al 2011), whereas the remaining objects in the local 
field have only AIS data. One more source of uncertainty in 
determining the SFR using SED fitting is connected with the 
large errors introduced by the SDSS pipeline on the photometry 
of large nearby galaxies due to "shredding" in multiple pieces, 
which leads to wrong magnitude determinations. 
Given these considerations, we decided to make a new estimate 
of the "instantaneous" SFR based on our own measurements of 
the hydrogen recombination lines using the data from the Ha3 
narrow -band imaging survey of nearby galaxies selected from 
ALFALFA in the spring sky of the Local Supercluster (Gavazzi 
et al. 2012a, Paper I), including the Virgo cluster. 

Virgo is the cluster of galaxies nearest to us, which, in spite 
of the enormous angular extent of about 6x10 square degrees, 
lies in the footprint of both ALFALFA and Ha3. Already in 
early studies Virgo appeared to be a cluster in the early stage 
of dynamical evolution. Its general appearance is highly irregu- 
lar (Abell et al. 1989), and it is classified as Bautz-Morgan type 
III (Bautz & Morgan 1970). Owing to the intensive optical sur- 
vey by Binggeli et al. (1985), its velocity dispersion was found 
to be 760 + 45 km s _1 , lower than that of evolved clusters, in- 
dicating a shallow gravitational potential. The spiral fraction is 
48% compared to 16% of Coma. As seen by ROSAT (Bohringer 
et al. 1994), the Virgo cluster is known to have an irregular X- 
ray morphology. These features provide evidence that the Virgo 
cluster is currently undergoing dynamical evolution (Forman & 
Jones 1982). 

The role of the environment in shaping the star formation 
properties of galaxies in the local Universe is still an open 
question. As early as in the Las Campanas redshift survey 
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(Hashimoto et al. 1998), or in the DR3 release of the SDSS 
(Gomez et al. 2003), evidence has built that suppression of the 
SFR of galaxies in dense environments takes place at z=0, up to 
z=0.8 (Patel et al. 2009) (see a review by Boselli & Gavazzi 
2006). By comparing the statistical HI properties of galaxies 
in nine nearby clusters with those of "field" objects of similar 
type and luminosity, Giovanelli & Haynes (1985) showed that 
clusters contain a high percentage of HI deficient galaxies. The 
mechanism that most likely contributes to significant gas deple- 
tion in clusters is ram-pressure stripping (Gunn & Gott 1972). 
Since then the HI deficiency parameter has been used as a proxy 
for the degree of the perturbation that galaxies experience in the 
harsh environment of rich galaxy clusters. 

Taking advantage of the presence of the Virgo cluster in Ha3, 
the first goal of the present work is to quantify, in addition to the 
HI content, how the SFR is affected by the environment of an 
evolving cluster such as Virgo, following the line traced by the 
pioneering work of Kennicutt (1983). We approach this question 
by studying a possible residual correlation between the SFR, the 
gas content and the gas deficiency after the first-order scaling 
law from the stellar mass have been removed. 

The dataset used in this paper has been presented in Paper I 
of the present series, which illustrates Ha3, the ongoing imag- 
ing survey at the 2.1m telescope of the San Pedro Martir (SPM) 
Observatory, the details of the sample selection, completeness, 
the Ha observations, data reduction and analysis, which gives 
access to the integrated extinction-free measurements of the SFR 
for 233 galaxies observed at SPM between 2006 and 2009. The 
present Paper II contains the analysis of the integrated SFR ob- 
tained using Ha3. After a brief introduction to the galaxy sam- 
ple (Section |2]), we derive an updated calibration of the HI defi- 
ciency parameter valid for dwarf gas-rich systems such as those 
that dominate Ha3 (Section |3j. The main results of the present 
analysis are given in Section [4] where we analyze the relations 
between the mean optical colors and the extent of the star for- 
mation regions of galaxies in the present survey as a function of 
their HI content. The scaling relations between gas and stellar 
masses and current SFR in the Local Supercluster are analyzed 
(Section|5j- The overall perturbation on the gas content and star 
formation properties of galaxies in the Local Supercluster due 
to the Virgo cluster is quantified in Section [6] The results are 
discussed in Section [7] where some conclusions on the present 
evolutionary stage of the Virgo cluster are drawn, and we sum- 
marize our results in Section [8] 

Two companion papers continue this series. Paper III (Gavazzi 
et al. 2012b) carries an analysis similar to the one reported in the 
present paper, focused on the properties of galaxies in the Great 
Wall, including the Coma cluster. Paper IV (Fossati et al. 2012) 
contains the analysis of the structural properties of galaxies in 
the local and Coma Superclusters. The CAS parameters (con- 
centration, asymmetry, and dumpiness developed by Conselice 
2003) are determined on both the r band and the Ha images to 
study the morphology of the star-forming regions and compare 
it with that of older stars at arcsec resolution (unachievable on 
GALEX UV data which have a resolution of 4-5 arcsec). 

Throughout the paper we adopt a flat ACDM cosmology, 
with Ho = 73 km s _1 Mpc~' and Q.^ = 0.7. All magnitudes are 
given in the AB system unless explicitly noted. 

2. Sample 

The sample analyzed in this work is drawn from the 900 square 
degree region ll h < R.A. < 16 h ; 4° < Dec. < 16°, covering 
the Local Supercluster, including the Virgo cluster. This region 



has been fully mapped by ALFALFA which provides us with a 
complete sample of Hi-selected galaxies with gas masses as low 
as 10 7 - 7 Mq (Haynes et al. 201 1). 

Ha3 is a follow-up survey consisting of Ha imaging obser- 
vations of ALFALFA targets with high signal-to-noise (typically 
S/N > 6.5) and a good match between two independent po- 
larizations (code = 1 sources; Giovanelli et al. 2005, Haynes 
et al. 2011). Ha observations of members of the Virgo cluster 
were given in Gavazzi et al. (2002a, b, 2006), Boselli & Gavazzi 
(2002), and Boselli et al. (2002b). Images and fluxes are also 
publicly available via the GOLDMine web server (Gavazzi et 
al. 2003). Ha data for galaxies in low-density regions of the 
Local Supercluster are given in Paper I. Outside Virgo the red- 
shift window of Ho-3 is 350 < cz < 2000 km s _1 , and we lim- 
ited the study to objects with HI fluxes Fm > 0.7 Jykms -1 , 
while in the Virgo cluster the velocity interval is extended to 
350 < cz < 3000 km s to map the cluster in its full extent 
we did not restrict ourselves to Fm > 0.7 Jy km s . At the dis- 
tance of 17 Mpc, assumed for the Virgo subcluster A (Gavazzi 
et al. 1999), a flux limit S21 = 0.7 Jy km s -1 corresponds to an 
HI mass M ffl = 10 7 7 M Q , while Ha3 is complete to M m > 10 8 
M Q (see Paper I). Figure [T] illustrates the sky region covered by 
Ha3, which contains 409 galaxies. 

The basic quantities extensively used in this paper, i.e, the HI 
mass, the stellar mass, and the extinction-free SFR were derived 
in Paper I. The HI mass is defined by Mhi = 2.36 ■ 10 5 • S21 ■ 
D 2 , where D is the distance to the source in Mpc and S 21 is 
the integrated flux of the HI profile in units of Jy km s _1 from 
ALFALFA. 

As explained in Paper I, the g and i band total magnitudes of 
galaxies detected by ALFALFA were measured by us on SDSS 
images (DR7, Abazajian et al. 2009) using IRAF. This ensured 
proper foreground star-subtraction and overcomes the "shred- 
ding" problem generated by the automatic SDSS pipeline. 

Inside the Virgo cluster we obtained an optically selected 
sample (see SectionUJi of 640 objects based on the VCC catalog 
by Binggeli et al. (1985) (including objects of all morphologi- 
cal types with B < 18 mag between 4" < Dec < 16° and ex- 
cluding those detected by ALFALFA and galaxies in the back- 
ground of Virgo with cz >3000 km s~'). For this large sample 
we computed the integrated magnitudes using an automatic pro- 
cedure. Starting from the SDSS images that are available from 
the GOLDMine web site, we produced deep "white" images 
by stacking frames in five filters. Using SExtractor (Bertin & 
Arnouts 1996), we separated the galaxies of interest from fore- 
ground stars and background galaxies. After masking these ob- 
jects, we computed the integrated magnitudes using elliptical 
apertures, the parameters of which were computed on the white 
images. As a test, we compared the photometry derived with this 
automatic procedure with a subset of galaxies that were mea- 
sured manually on i and g images (Paper I) and we obtained 
very consistent results. 

g - i colors were corrected for internal extinction using the 
empirical transformation dependent on stellar mass and galaxy 
inclination, derived in Appendix A of Paper III: 

(g-0 o = (g-0-{+0.17-[l-cos(md)]-[log(^)-8.19]}, (1) 

where incl is the galaxy inclination (computed for disk galaxies 
(from Sa to Sdm) following Solanes at al. (1996). 

The stellar mass was derived from the i magnitudes and (g - 
i) color using the transformation 



log© = -1.94 + 0.59- (,- 0o + 1.15- logg) 



(2) 
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Fig. 2. Sky distribution of galaxies around the Virgo cluster. ETGs (from the VCC with m p <18.0 mag) are given with small black 
dots. Hi-selected LTG galaxies, coded according to their HI deficiency parameter (red dots: Defm > 0.9; green dots: 0.3 < Defm < 
0.9; blue dots Defm ^ 0.3). The empty red circles represent the 26 LTGs undetected by ALFALFA listed in Table Q] Superposed 
are the X-ray contours from ROSAT (Bohringer et al. 1994). The position of cluster A (M87), B (M49) and of clouds M, N and W 
are given. Clouds M and N coincide in position but have significantly different distance moduli. 



consistent with the mass determination of MPA-JHL0, where 
log Li is the i band luminosity in solar units (log = (7-4.56)/- 
2.5). 

The extinction-free SFR was computed from the luminosity 
of the Ho- line after correcting for Galactic and internal dust ex- 
tinction and deblending from [Nil] emission. Details on these 
quantities are given in Paper I. 

Following a similar procedure as in Gavazzi et al. (2010, 
2011) and Paper III, the local surface density S around each 
galaxy was computed within a circle of 1 Mpc radius. Around 
each galaxy we computed the density contrast as 

I- < E > 



< S > 



, where S is the local surface density and < S > = 1.3 gal 
(Mpc) -2 represents the mean surface density measured in the 
whole region. 

3. HI content of Virgo galaxies 

The Virgo cluster is well known for its unrelaxed structure 
(Forman & Jones 1982) and for its irregular morphology, elon- 
gated in the N-S direction. Figure [2] shows an enlargement of 
Figure [TJ highlighting the zone that contains the Virgo clus- 
ter, with the early-type galaxies (hereafter ETGs) added to em- 
phasize the global density of galaxies in this region. Following 



n 






Def„, 



Fig. 3. Histogram of the Defyi parameter separately for the iso- 
lated and W, N clouds (black) and for members of the Virgo A 
and B subclusters (red). The bulk of Hi-deficient galaxies is lo- 
cated within the two subclusters. 



www.mpa-garching.mpg.de/SDSS/DR7/, see Salim et al. (2007) 
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Gavazzi et al. (1999), we assume that Virgo is composed of the 
main cluster A, which contains the bulk of the X-ray emission 
centered on M87, and a secondary cluster B centered on M49, 
infalling into A at about 750 km s _1 . These subcomponents are 
at a distance of about 17 Mpc. Clouds M and W are twice more 
distant and are probably infalling into the main cluster along 
with cloud N, which is in the process of merging with cluster 
A. Virgo is a spiral -rich cluster. The Virgo Cluster Catalogue 
(VCC, Binggeli et al. 1985) lists 48% (36%) galaxies of type 
later than SOa for m p < 18 mag (m p < 20). Restricting this to 
trip < 13.9 mag, corresponding to the completeness limit of the 
spectroscopic SDSS g < 17.75 at Coma, the spiral fraction in- 
creases to 60%, i.e, four times more than in the Coma cluster 
(16%) at similar depth. 

It might well be that some of the late-type galaxies (hereafter 
LTGs) that contribute to the large spiral fraction in Virgo are 
LTGs seen projected on Virgo, but that are distributed along the 
line of sight on a filament connecting the local group to Virgo, 
similar to some of the elongated filaments that exist in the Coma 
Supercluster (see Paper III). However, we estimate that if we 
were to look at Coma from an angle such that we integrate along 
one of these filaments, the number of LTGs projected on Coma 
would increase from 16 % to about 30%. Although significant, 
this is not enough to explain the high fraction of LTGs seen in 
Virgo in terms of projection effects. 

The parameter that is bestsuited to disentangle galaxies that 
physically belong to the Virgo cluster from either projected on 
it or belonging to clouds not yet fully processed by the cluster 
potential is perhaps the Hi-deficiency parameter (Defni). High- 
resolution HI maps of Virgo galaxies, such as those obtained by 
the VIVA survey of Chung et al. (2009), confirm earlier claims 
(Cayatte et al. 1990) that the perturbed spiral/Irr galaxies in the 
Virgo cluster have a high Defm and their HI morphology shows 
evidence of truncation/one-sidedness, consistent with the hy- 
pothesis that they are suffering from ram-pressure stripping by 
the dense cluster ICM. 

The Defm parameter has been defined by Haynes & 
Giovanelli (1984) as the logarithmic difference between the HI 
mass observed in a galaxy and the expected value in isolated and 
unperturbed objects of similar morphological type T and linear 
diameter d^. Defm =< log Mhi(T, df) > - logMniiobs). Here, 
< log Mhi(J, di) >= C\ + C2 x 2 log(d[}, where d^ (in kpc) is 
determined in the g band at the 25' /! mag arcsec~ 2 isophote. The 
coefficients C\ and C2 were determined by Haynes & Giovanelli 
(1984) by studying a control sample of isolated objects, and 
later, on a larger sample by Solanes et al. (1996). Both sam- 
ples are composed almost exclusively of giant spirals however. 
The Hi-deficiency parameter was therefore poorly calibrated for 
dwarf objects. 

The problem is re-addressed here using ALFALFA by as- 
suming as isolated the galaxies outside the Virgo cluster (out- 
side the interval 12 /! 08'W < RA < 12 ft 48 m 00 s ), excluding a 
few more objects that appear to be clustered in small groups. 
We took these objects (the blue symbols in Figure [TJ as a ref- 
erence sample of nearly isolated objects, representative of local 
normal and unperturbed galaxies. Scd-BCD dwarf systems dom- 
inate this sub-sample. 

The coefficients C\ and C2 obtained by us for the giant galax- 
ies (Sa-Sc) are similar to those of Solanes (1996) and of Boselli 
& Gavazzi (2009). Consistent values are found in our sample 
for Scd-BCD. For this reason, we decided to adopt C\ = 7.51 
and C2 = 0.68 for all LTGs (Sa-BCD), obtained by combining 
all late-type galaxies. The resulting Defm are plotted in Figure 
[3] where we display the frequency distribution of Defm in the 



isolated sample mixed with members of Virgo clouds M and 
W (black histogram), separately from Virgo members of sub- 
clusters A and B (red). For the isolated and M, W clouds, we 
find a mean deficiency < Defm >= 0.01 + 0.35, consistent 
with normal HI content, while for sub-clusters A and B this is 
< Defm >= 0.33 + 0.39. 

Unfortunately, the two distributions have overlapping tails, 
which reflecting the difficulty of spatially separating in 3-D) 
galaxies that belong to the real cluster from those that are just 
projected on it. For example (see Figure |2} there are many per- 
fectly Hi-normal LTGs projected within 2 degrees from M87, 
while there are several highly deficient objects up to 4 degrees 
away from M87 or outside the X-ray contours. However, the 
large majority of LTGs outside the X-ray contours are non- 
deficient (blue symbols). Thus the assumption that the Hl- 
deficiency parameter is a reliable tracer of the degree of pertur- 
bation suffered by LTG by the cluster medium, allows us to adopt 
in the rest of this paper a conservative threshold of Defm=03 
(1 cr above the mean deficiency of the normal sample) to sepa- 
rate normal from deficient galaxies. There are approximately 30 
galaxies in the bin 0.2 < Defm < 0.3 that we consider normal, 
because they are projected onto the N and W cloud. In reality 
they might belong to the outskirts of the A cluster, which would 
justify their slightly positive deficiency. Throughout this paper 
we adopt Defm = 0.3 as a threshold below which (270) galax- 
ies can be considered as normal or unperturbed by the cluster 
environment, while those with Defm > 0.3 are treated as envi- 
ronmentally perturbed. The abundant number statistics allows us 
to further sub-divide the deficient sample into 0.3 < Defm < 0.9 
(115 perturbed), and Defm > 0.9 (24 highly perturbed). 

4. Radio versus optical selection 

The comparison between the star formation properties of galax- 
ies selected with radio vs optical criteria helps to shad some light 
on the relation between the star formation properties (both global 
and nuclear) of LTGs and their HI content in various environ- 
ments. 

As already stressed in Paper I, the radio selection biases Ha3 
toward gas-rich, late-type galaxies (see also Gavazzi et al. 2008). 
The color (g-i)„ versus stellar mass diagram composed of galax- 
ies in Ha3 (top panel of Figure |4]i contains no well-developed 
red sequence (as emphasized by the red dashed line that marks 
the loci potentially occupied by the red sequence). This disagrees 
with similar diagrams derived from optically selected samples 
(e.g. Hogg et al. 2004). The census of Ha3 is made almost ex- 
clusively of LTG galaxies in the blue sequence. Only 15 ETG 
galaxies are detected in the radio in our survey (see Figure top 
panel). Interestingly, there is a minor (0.13 mag) but significant 
difference in the mean {g - i) a color between LTGs with normal 
HI content (blue symbols in Figure [4] top panel) and LTGs with 
deficient HI content (cyan symbols in Figure|4l top panel). 
Furthermore, not all LTGs are detected by ALFALFA. To 
demonstrate this point one would like to be able to extract 
a priori an optical catalog of galaxies in the sky covered by 
ALFALFA. A similar exercise was carried out by Cortese et al. 
(2008a) in their HI survey of the cluster A1367 and in Paper III 
for the Coma cluster. 

Owing to the notorious limitations that affects the completeness 
of SDSS for z « 0.05 due to the shredding problem (Blanton et 
al. 2005), a similar exercise cannot easily be repeated for the 
Local Supercluster and Virgo. Alternatively, we proceeded in 
two steps as follows. Out of the boundaries of the Virgo clus- 
ter we searched in NED all galaxies with cz < 2000 km s _1 
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in the footprint of Ha-3. At their position we searched the SDSS 
navigator and we took the i and g Petrosian magnitudes. After se- 
lecting those with i <17.0 mag, we were left with 50 objects, 35 
of late and 15 of early type. It is not surprising that the 15 ETGs 
are not detected by ALFALFA. The 35 LTG have stellar mass 
logM* < 8.5, thus it is not surprising that they remain undetected 
by ALFALFA because, owing to the relation logMni = 0.55 x 
logM, + 3.85 that exists for unperturbed galaxies (see Section [8] 
and Gavazzi et al. 2008), they lie below the survey limiting sen- 
sitivity. 

Inside the Virgo cluster a reliable optical selection was obtained 
using the VCC catalog by Binggeli et al. (1985) (including ob- 
jects of all morphological type with B < 18 mag (Vega) between 
4° < Dec < 16° and excluding those detected by ALFALFA and 
galaxies in the background of Virgo with cz >3000 km s~'). 

For this large-sample we computed the integrated magni- 
tudes using the automatic procedure described in Section |2] In 
the surveyed region we count 640 optically selected galaxies 
with r <17.5 mag (in addition to 509 HI selected galaxies). The 
color-mass relation of optical selected galaxies is plotted in the 
bottom panel of Figure [4] separately for early (red symbols) and 
late types (green symbols). The 26 massive LTG Virgo members 
(log(M»/M ) > 9) undetected LTGs even though they are in 
the footprint of ALFALFA0 are listed in Table [1] They are well- 
known high-HI-deficiency objects from previous deep pointed 
observations taken at Arecibo by Helou et al. (1984), Haynes & 
Giovanelli (1986), Hoffmann et al. (1987, 1989), and Gavazzi et 
al. (2005). The interesting new fact is that the undetected LTGs 
have redder colors than the corresponding Hi-detected galaxies 
(see bottom panel of Figure 0}, and are approximately 0.2 mag 
bluer than galaxies in the red sequence, which fill the "green 
valley" (see also Cortese & Hughes 2009). 

All fits have a consistent slope of 0. 10 and the intercepts are: 
0.18 for ETGs; -0.05 for undetected LTGs (0.23 mag bluer); - 
0.24 for detected LTGs with Def HI > 0.3 (0.19 mag bluer); and 
-0.37 (0.13 mag bluer) for detected LTGs with Def HI < 0.3. A 
consistent increase of the specific star formation rate (SSFR) in 
steps of decreasing Deffii is also visible in Figure [8] (c). The 
above exercise, although it should be taken with a grain of salt 
because of NED incompleteness, confirms the result of Cortese 
et al. (2008a) (A1367) and of Paper III (Coma): given the sensi- 
tivity of ALFALFA at the distance of the Local Supercluster, 
only LTGs are detected. Among them, however, ALFALFA 
misses objects that are the most Hi-anemic. 
Another way to show how the morphological/HI-content mix of 
galaxies depends on the environment in and around the Virgo 
cluster is offered in Figure [5] Here we plot the corrected color of 
galaxies of all morphological types regardless of their HI proper- 
ties (i.e, including objects undetected by ALFALFA) as a func- 
tion of the projected separation from M87 (bottom panel) and of 
the local galaxy density (top panel). We use the same symbols as 
in Figure|4] The Virgo cluster (which lies within 10 degrees from 
M87) is visible because it causes a strong color and morphology 
segregation, i.e, all ETGs (in the red sequence) are confined into 
its boundaries. Moreover, it also segregates LTGs with a low HI 
content from their HI -rich counterparts, which are mostly found 
outside its boundaries. These families also appear to be signif- 
icantly segregated as far as their specific star formation is con- 
cerned, as indicated by their different color. The properties of the 



3 Six additional VCC objects are undetected by ALFALFA, not be- 
cause the are Hl-poor, but because they are confused among each-other 
or with other sources: i.e, VCC 386, 483, 497, 1673, 1676, and 1972; 
see Hallenbeck et al. (2012). 



1.5 



OB 

a 



0.5 


1.5 



so 

CO 



0.5 



F-WrfT*. .. * LTG Def HI S0.3 

_^ # * . . • aa Det. 



J I M 



I I I I 



MM 



ETG 

I I I I I I I I I 



+4 




_L 



aa Undet. 



_L 



LTG 
ETG 



9 10 
logM. (M e ) 



1 1 



12 



Fig. 4. Top panel: (g - /)„ color (corrected for internal ex- 
tinction) vs. stellar mass diagram of Ha3 galaxies detected 
by ALFALFA (isolated-H Virgo), color-coded by morphology 
(red=ETGs, biue+cyan=LTGs). Bottom panel: color-mass dia- 
gram of galaxies not in Ha3 (undetected by ALFALFA). They 
have been optically selected from NED among isolated ob- 
jects, and from the VCC within the Virgo cluster (red = ETG; 
green=LTG). In both panels the diagonal dashed red lines mark 
the loci of the red sequence. The 26 LTGs in the VCC, unde- 
tected by ALFALFA, with logM„ > 9.0 are listed in Table Q] 



studied galaxies on nuclear / circumnuclear scales are investi- 
gated in Figure|6]and in Table|2] The figure gives the cumulative 
distribution of the ratio of total Ha flux (measured in the imag- 
ing material) to the Ho- flux measured in a central aperture of 
variable size according to the total size of the galaxy, separately 
for LTGs-in ALFALFA with Def m < 0.3 (blue) and for LTGs 
with Defyi > 0.9 mixed with LTGs undetected by ALFALFA 
(red). A Kolmogorov-Smirnov test indicates that the probability 
that the two samples are drawn from the same parent population 
is 0.03 %, i.e, that they significantly differ. In other words, the 
Ha sources associated with Hi-rich galaxies are extended, while 
those of strongly deficient/undetected objects are predominantly 
nuclear/circumnuclear, see also Koopmann & Kenney (2004a,b). 
Table [2] summarizes the classification of the available nuclear 
SDSS spectra of LTG galaxies in this work in three bins of Defm 
(the LTG galaxies undetected by ALFALFA are included in the 
highest deficiency bin). It shows that the fraction of star-forming 
nuclei decreases significantly with increasing Defuu although 
some signs of nuclear star formation remain in 46 % of the most 
Hl-poor systems. The fraction of passive spectra is relevant only 
in the highest Defm bin, while the fraction of AGNs (including 
LINERS) does not significantly change with Defui (consistently 
with Gavazzi et al. 201 1). 
Altogether we find evidence that 

- there is a mild gradual progression toward redder colors from 
(1) LTGs detected by ALFALFA with normal HI content, (2) 
LTGs detected by ALFALFA with deficient HI content, (3) LTGs 
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Fig. 5. Corrected (g - /)„ color of galaxies of all morphological 
types as a function of the local galaxy density (top panel) and 
of the projected angular separation from M87 (bottom panel). 
Hi-rich LTGs are depicted with blue symbols; Hi-poor LTGs are 
plotted with cyan symbols; LTGs undetected by ALFALFA are 
given with green symbols; ETGs are marked with red symbols. 



Table 1. 26 LTG VCC galaxies in the Virgo cluster with logM, > 
9.0 not detected by ALFALFA, with detections or more stringent 
upper limits from previous pointed observations. A classification 
of the nuclear spectrum is given in Column 6 and the morphol- 
ogy of the Ha source in Column 7. 
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Fig. 6. Cumulative distribution of the ratio of total Ha flux in the 
imaging material to the Ha flux measured in a central aperture of 
10, 15, 30 or 45 arcsec diameter for galaxies of major axis 0-100, 
100-200, 200-300 and > 300 arcsec. The blue histogram refers 
to LTG-in ALFALFA with De fni < 0.3 and the red one includes 
LTGs with Def H i > 0.9 and LTGs undetected by ALFALFA. 



undetected by ALFALFA and (4) ETGs; 

- the same morphology/HI-content sequence depends on the am- 
bient density; 

- the extent of the star formation region in Hi-rich (1) LTGs is 
statistically larger than in their Hl-poor (2+3) counterparts. 

In other words, we confirm the finding of Paper III (Coma) that 
the environmental conditions not only affect the HI content of 
galaxies, but also govern the average color (specific star forma- 
tion rate) of galaxies. Moreover, whatever ablation mechanism 
is causing this pattern, it must proceed outside-in, in agreement 
with Boselli et al. (2008). 



Table 2. Frequency of nuclear activity in LTGs detected by 
ALFALFA with logM, > 7.0 (blue symbols in Fig |4j divided in 
three classes of increasing Defm- The undetected LTGs (green 
symbols in Fig|4| are combined with the Defni >0.9. The last 
column gives the percentage of available nuclear spectra. 



Def HI 


% PASS 


% AGN+LIN 


% HII 


compl.% 
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0.3-0.9 


14 


13 


73 
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>0.9 


24 


15 
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78 



5. Scaling relations 

The study of the scaling relations among the barionic con- 
stituents of (gas rich) galaxies has made a significant step for- 
ward with the work of Huang et al. (2012), who analyzed a.40 
(Haynes et al. 201 1), the most updated ALFALFA catalog avail- 
able so far. Here we compare our local data with the results of 
Huang et al. (2012), whereas in Paper III we compare Virgo with 
Coma. 

We begin by comparing in Figure [7] the star formation esti- 
mates obtained in the present work from the Ha corrected lu- 
minosities taken from Paper I, with the estimates of Huang et 
al. (2012) computed by fitting UV, u,g,r,i,z SEDs. The corre- 
lation cefficient is 0.65, implying that the two quantities have a 
probability higher than 99 % of being correlated. The bottom 
panel shows the distribution of the differences of the two quan- 
tities. The Gaussian fitted to the data was obtained after sigma 
clipping of the most deviant points. Two most extremes outlier 
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galaxies (VCC 131 and VCC 1932) merit a more detailed com- 
ment. These are both bright edge-on galaxies with significant 
Ha flux (see GOLDMine) (logF=-12.68 and -12.32 erg cirrV 1 
respectively), therefore with consistent SFR as computed in this 
work (-1.15 and -0.37 M Q yr _I ). Furthermore the SFR derived 
by directly converting the FUV magnitudes (corrected follow- 
ing Cortese et al. 2006, 2008b, assuming balance between the 
absorbed UV photons and the total FIR energy) into SFR are 
found to be (-0.44 and 0.3 1) within a factor of 5 from each other. 
Conversely, the SFR found by Huang et al. (2012) (1.11, -4.42) 
differs by more than five orders of magnitude, revealing some 
problems that are probably connected to the SED fitting proce- 
dure. Similarly, the galaxy VCC 785 is discrepant in Figure |7] 
the SFR by SED fitting is approximately two orders of magni- 
tude above the linear correlation. This is a face-on galaxy sur- 
rounded by a faint outer disk (as revealed in a deep 1500 sec 
NUV image by GALEX not accompanied by a similar length 
FUV exposure). However, as explicitly mentioned by Huang et 
al. (2012), when deep GALEX images are not available in both 
NUV and FUV bands, the single long exposure is discarted and 
the SED fitting is performed with available FUV and NUV data 
of similar exposures, albeit shorter. For VCC 785 the GALEX 
database reports a FUV magnitude of 19.5 mag compared with 
a NUV of 15.4 mag (as derived from 200 sec exposures). The 
SED fitting routine probably assumes an overestimated extinc- 
tion correction to allow for such a large discrepancy in magni- 
tude, hence it overestimates the SFR. In summary, we conclude 
that, albeit the large scatter, the two sets of SFR measurements 
are consistent, which makes the comparison between the scaling 
relations obtained from the two SFR estimators meaningful. 

The scaling relations relevant to our study are given in Figure 
[8] (a) through (c), and the linear regression parameters are listed 
in Table Hf. 

We consider the relation between the HI and stellar mass by 
showing in Figure[8](a) the HI gas fraction Mm /M* as a function 
of the stellar mass. The diagonal dot-dashed line represents the 
sensitivity limit of ALFALFA, computed for galaxies with an in- 
clination of 45 deg in the plane of the sky0. By slicing the Ha3 
sample in three subsamples of increasing Defni'. Defni < 0.3 
(normal; blue symbols); 0.3 < Defni < 0.9 (moderate defi- 
ciency; green symbols), and Defni > 0.9 (high deficiency; red 
symbols) we explore the whole range of environmental perturba- 
tions exhibited by galaxies in the Local Supercluster, an aspect 
that was not examined by Huang et al (2012). 
Consistently with the results of Gavazzi et al. (2008) and Huang 
et al (2012), the gas to stars fraction decreases by approxi- 
mately 4 orders of magnitude by increasing galaxy mass, from 
logM, ~ 7 M to ~ 11.5 M , independently of the HI de- 
ficiency parameter. This basic result seems to be independent 
of the fact that Hor3 is an Hl-selected survey. Indeed our scal- 
ing relation perfectly agrees with the one found in the GASS 
survey by Catinella et al. (2010) (see the black symbols in 
Figure |8ja)). Similarly, over a broader interval of stellar mass 
10 8 - 5 < M* < 10 11 5 M there is also a good agreement with the 
relation found by Cortese et al. (201 1) among Hi-normal galax- 

4 Linear regressions are obtained in this work using the "bisector" 
method by Isobe et al. (1990) (mean coefficients of the direct and the 
inverse relation). 

3 As introduced in Giovanelli et al. (2005), ALFALFA is a noise- 
limited survey rather than a flux-limited one. At any given integrated HI 
mass the 21 cm flux per velocity channel is inversely proportional to the 
width of the HI profile, hence to the galaxy inclination. The complete- 
ness and sensitivity of ALFALFA are clearly defined and discussed in 
detail in Saintonge (2007), Martin at al. (2010) and Haynes et al. (201 1). 
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Fig. 7. (Top) Correlation between the SFR used in this work 
(from Ha) and those adopted by Huan et al. (2012) (computed 
from SED fitting) The one-to-one relation is given with a dot- 
ted line. Some outliers marked with their VCC names are com- 
mented on in the text. (Bottom) The distribution of the residual 
of the two star formation indicators, with a Gaussian of mean fi 
and cr fitted to it. 



ies in a near-infrared selected sample derived from the Herschel 
Reference Survey (HRS; Boselli et al. 2010) (see the large ma- 
genta squares in Figure |8la)fl This result is consistent with 
the downsizing scenario (Gavazzi 1993, Gavazzi et al. 1996, 
Gavazzi & Scodeggio 1996, Boselli et al. 2001, Fontanot et 
al. 2009), where progressively more massive galaxies have less 
gas/star at z=0 because they have transformed most of their gas 
into stars at higher z, while late-type dwarf galaxies retain large 
quantities of hydrogen capable of sustaining the star formation 
at some significant rate at the present cosmological epoch. 

The logSFR versus logM, relation (Figure [8jb)) obtained 
at z=0 is near the direct proportionality (slope=0.97). A nearly 
exponential mass growth has also been obtained at redshift up 
to 3 by, e.g., Griitzbauch et al. (2011) and beyond, e.g. Labbe 
et al. (2010). The specific star formation rate (SFR/M,) (Figure 
|8jc)) is thus unsurprisingly flat, or marginally decreasing with 
increasing mass. 

In none of Figures [8] there is compelling evidence that the 
relations shown have a significant change of slope in the con- 



6 The GASS and HRS points plotted in Figure (8{ a) are binned aver- 
ages, thus the error-bars do not represent the dispersion of the distribu- 
tion. 
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sidered interval of stellar mass, in contrast to the one found by 
Huang et al. (2012) near logM* ~ 9 M . This is probably be- 
cause the Local Supercluster lacks high mass galaxies due to the 
small sampled volume compared to a AO. However, a change of 
slope 

Good agreement is found between Ha3 and HRS (Cortese 
et al. 201 1) by comparing our deficient (green dots) or strongly 
deficient (red dots) objects with the Virgo HRS subsample (large 
red squares in Figure [8j a)). Altogether, there is full agreement 
with previous surveys that on the common mass range there is 
a primary scaling relation between HI and the stellar mass, but 
Ha3 extends the studied range of mass to one order of magnitude 
fainter. 

High-deficiency galaxies (see Figure |8la)) are found exclu- 
sively among the most massive systems and lie dangerously near 
to the ALFALFA detection limit already at the distance of Virgo 
(this lack of sensitivity is at the origin of the LTGs undetected by 
ALFALFA considered in Section More highly Hi-deficient 
galaxies exist at low stellar mass, but they require deeper HI 
observations than ALFALFA to be detected. Indeed when the 
optically selected galaxies (see Section |4]l are considered and 
their deeper available HI observations are used to complement 
the ALFALFA data in Figure [HJa), several objects with severe 
HI depletion are detected at all stellar masses, including objects 
as faint as 10 7 M Q . This is consistent with the idea that hydrogen 
depletion occurs more easily in dwarf galaxies within a cluster 
like Virgo (Boselli et al. 2008). The time scale for gas deple- 
tion is short enough (~ 100 Myr) to quickly transform them into 
anemics that become gas-free, passive dEs. This is precisely the 
mechanism invoked by Boselli et al. (2008) and by Gavazzi et 
al. (2010) to migrate dwarf star-forming systems into dEs in the 
neighborhood of rich galaxy clusters. 



Table 3. Linear regression coefficients of the relations in Figures 
[8] obtained using the bisector method (mean coefficients of the 
direct and the inverse relation). 
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In our sample the specific star formation is found to depend 
only mildly on the stellar mass (Figure |8jd), compared with 
Huang et al (2012) or with optical selected samples (Gavazzi 
et al. 2002b, Schiminovich et al. 2010). 

The most remarkable feature of the specific star formation 
vs stellar mass relation is the enormous scatter, as noted by other 
authors (Lee at al. 2007, Bothwell et al. 2009), particularly ev- 
ident at low mass, where the scatter reaches 2 orders of mag- 
nitude. The faintest dwarfs are amorphous objects composed of 
few HII regions. Some of these are outliers with respect to the 
trends in Figures|8] with a high SFR for a given HI mass and M». 



It is debated (Boselli et al. 2009, Fumagalli et al. 201 1, Weisz et 
al. 2011) whether sporadic bursts of SFR occurring in the indi- 
vidual HII regions in a stochastic manner could produce 100 % 
variations in the integrated SSFR, hence a spread similar to the 
one observed in Figure |8fc). The same relation was also stud- 
ied by Bothwell et al. (2009) and by Schiminovich et al. (2010) 
(see their Fig. 8) in the (optically selected) GASS survey, lim- 
ited to galaxies with 10 10 < M, < 10 1L5 M . In this mass range 
our data agree with their high SFE objects (see the blue line 
in Figure [He)), but in addition to these, Bothwell et al. (2009) 
and Schiminovich et al. (2010) have many low SFE objects at 
log SSFR ~ -12 yr~'(see the red line in Figure|8lc)) that are sim- 
ply not sampled in Ha3, which makes the overall relation steeper 
than in Ha3. Boselli et al. (2001) also found a correlation steeper 
than in Figure Oc). 

Summarizing, our analysis has highlighted in agreement 
with previous studies (e.g., Vulcani et al. 2010, Sobral et al. 
2011) that the stellar mass is a fundamental parameter in driv- 
ing the observed scaling relations. The effect of the environment 
is then superimposed on this scaling relations and is more evi- 
dent in the relations that involve the mass of neutral hydrogen, 
which more easily prone to the effects of the environment when 
quantified using the HI deficiency. 

6. Overview of the SFR in the Local Volume 

We now illustrate the global role of the Virgo cluster on the star 
formation properties of the Local Supercluster as a whole, as 
shown in Figure [9] Only the population of LTGs is considered, 
meaning we disregard ETGs that would obviously trace the phe- 
nomenon of morphology segregation in the Virgo cluster. The 
top panel shows the celestial distribution of LTGs in the present 
work. Annuli centered on M87 are drawn with dashed lines, cho- 
sen to contain approximately an equal number of galaxies^. After 
discarding galaxies that do not properly belong to cluster A or 
B (background groups, such as clouds M, W), within each an- 
nulus we compute the median specific SFR (SFR per unit stellar 
mass) (second panel). The depletion or Robert's time (HI mass 
divided by SFR) shown in the third panel gives the duration of 
the remaining star formation activity at the present rate for the 
available HI content. The inverse of the latter quantity is some- 
time called the star formation efficiency. The HI deficiency is 
given in the bottom panels. These quantities are plotted in bins 
of projected distance from M87 (top panels) and in bins of local 
galaxy density (bottom panels). 

There is a significant radial dependence of the HI content 
(the strongest of all), of the specific star formation rate, and 
therefore of the depletion time. LTGs in the Local Supercluster 
located up to 20 degrees away from Virgo have a life expecta- 
tion longer than 10 10 yr (i.e, almost one Hubble time) if they 
continue to burn their hydrogen at the present rate, while galax- 
ies near to M87 will remain active 10 9 3 yr (less than 2 Gyr). This 
derives from the reduced gas reservoir of galaxies in the cluster 
environment (Hi-deficiency) and appears to be associated with a 
significant suppression of the star formation rate. 

7. Discussion 

The scaling relations between the gaseous (HI) mass and the star 
formation rate as a function of the stellar mass, in intervals of HI 

7 Figure [9] was obtained adding to the ALFALFA-detected LTGs the 
undetected LTGs that were optically selected with HI and Ha measure- 
ments taken from GOLDMine. 
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Fig. 8. Scaling relations of HI mass, stellar mass and SFR for galaxies in three classes of increasing Defui- blue: Defni < 0.3 
(normal); green: 0.3 < Defni < 0.9 (perturbed), red: Defni > 0.9 (highly perturbed). Symbols are assigned according to the 
morphology. E-SO-SOa: circles; giant spirals (Sa-Sm): squares; Irr-BCD: triangles; AGN: asterisks. Panel (a): gas fraction (M////M*) 
versus M* . The diagonal dot-dashed line represents the limiting sensitivity of ALFALFA at the distance of Virgo. The black symbols 
give the representative - sample extracted from GASS by Catinella et al. (2010), the large magenta squares represent the Hi-normal 
sample and the large red squares the Virgo sample from Cortese et al. (201 1). Across the stellar mass range in common there is an 
almost perfect agreement, but the present sample extends the correlation to two orders of magnitude fainter. Panel (b): SFR versus 
M,. Panel (c): Specific SFR versus M„. The blue and red solid lines reproduce the relations plotted in Figure 8 of Schiminovich et 
al. (2010) for the star-forming sequence (blue) and the non-star forming locus of red galaxies (red). 



deficiency are analyzed in Figure [8] and are found to be consis- 
tent with previous determinations based on HI and on optically 
selected samples. 

The environmental conditions occurring in a young, still 
evolving cluster of galaxies such as the Virgo cluster have 
profound consequences on the morphological composition of 
its member galaxies, producing the well known phenomenon 
named morphology segregation (Dressel 1980). Moreover, it 
induces observable modifications to the population of LTGs, 
even when they are considered alone: their HI content and to 
a lesser extent their specific star formation are significantly re- 
duced compared to LTGs found in relative isolation. To quantify 
the quenching introduced by the environment on the gas con- 
tent and on the SFR independently of the stellar mass, we plot 
in the central panel of Figure [10] (a,b) the relation of the spe- 
cific star formation versus gas fraction. This Figure was obtained 
by dividing both the SFR and the HI mass by the stellar mass, 
so that both axes contain normalized quantities. In other words, 
this figure contains the residual correlation excluding the scaling 
law with stellar mass. The symbols are color-coded according to 
the De fiii parameter. The histograms in the side panels enable 
to estimate the average gas fraction and SSFR in three bins of 
HI deficiency. While the left panel contains only giant galaxies 
(logM* > 9 M ), the right panel shows all galaxies. Both pan- 
els show that even after subtracting the scaling law with stellar 
mass , a significant residual (second-order) dependence of the 
SSFR on the Hi-deficiency remains: gas-rich objects that con- 
tain as much mass in HI and stars (log Mhi/M„-Q) have SSFR 
ten times higher than galaxies that contain 100 times less HI 
relative to stars. The gas-rich, high SSFR systems are exclu- 
sively found among unperturbed galaxies (Defni < 0.3), while 
gas-poor objects with a lower SSFR are exclusively highly per- 
turbed (Defm > 0.9) members of Virgo. The comparison be- 
tween the two panels emphasizes the large scatter in the SSFR 
of low mass objects mentioned earlier. These faint systems pref- 
erentially have a high gas fraction, normal or intermediate HI- 
deficiency and an erratic specific star formation rate over 2 dex, 



in contrast to high stellar mass systems, whose specific star for- 
mation correlates better with the gas fraction (see Boselli et al. 
2009). In summary, once the scaling relation with stellar mass 
is neutralized, it is found that by decreasing the gas fraction by 
2 dex, the specific star formation rate decreases by one dex, i.e, 
the effect of the environment is ten times as effective at removing 
the HI gas than at quenching the star formation. 

This is consistent with the finding of Figure [9] which shows 
that the star formation efficiency (M#//SFR) decreases going to- 
ward the Virgo cluster, i.e, Mhi decrease more than SFR, im- 
plying that the cluster produces a significant perturbation to the 
Local Supercluster, even on its member late-type galaxies, par- 
ticularly efficient at reducing their HI content and to a lesser ex- 
tent at quenching their present star formation rate. 

In Section |4] we showed that there is evidence of a continu- 
ous sequence of transformations that make LTGs gradually gas- 
poorer and redder in their disk components while they approach 
the Virgo cluster. Meanwhile, their star formation shrinks to the 
circumnuclear or even nuclear scale. The last feature calls for a 
truncation process of the star formation that is inconsistent with 
the so-called strangulation (Boselli et al. 2006, 2008) and con- 
sistent with ram pressure. 

HI ablation occurs on very brief time scales of 100-200 Myr 
for low-mass galaxies and 300 Myr for massive galaxies. An 
example is NGC 4569, analyzed by Boselli et al. (2006). In spite 
of being very Hi-poor (85% of the HI gas has been removed from 
this galaxy by ram-pressure), neverless it is only marginally H2 
deficient (Fumagalli et al. 2009) and its global star formation is 
only reduced by a factor of two. This is consistent with the fact 
that most of the gas has been removed from the outer disk, while 
in the inner regions the star formation continues sustained by 
copious H2, coexistent with an inner HI shell. The time-scales 
for ram-pressure stripping in this massive galaxy is consistently 
estimated to be ~ 300 Myr by Vollmer et al. (2004) and between 
100 and 400 Myr by Boselli et al. (2006). Estimates of 200 Myr 
were derived for two dwarf galaxies in the Virgo cluster: VCC 
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Fig. 9. Top panel: the sky distribution of LTGs within seven an- 
nuli of increasing radius from M87. Galaxies in each ring are 
given with a different color, (b): SSFR; (c): SFR per unit HI 
mass; (d) Hi-deficiency as function of the projected angular sep- 
aration from M87. Bottom panel: same as above in bins of lo- 
cal galaxy density. This analysis is limited to galaxies with pro- 
jected distance < 45° from M87, i.e, 20 objects in the interval 
15:00:00< RA < 16:00:00) were excluded because they do not 
constitute a statistically meaningful subset. 



1217 (Fumagalli et al. 201 lb) and VCC 1249 in the M49 group 
(Arrigoni Battaia et al. 2012). 

The ongoing infall rate of galaxies on the Virgo cluster can 
be obtained by dividing the number of HI rich LTGs that still 
exist in the Virgo cluster by the time scale of the HI stripping 
mechanism. Table [4] contains two estimates: for > 10 7 M and 



for > 10 9 M Q (the second interval is to compare our estimate 
with a similar estimate obtained for the Coma cluster in Paper 
III). In the two cases we obtain a current infall rate of 400 - 100 
galaxy Gyr~' in the two mass bins (the former consistent with 
300 galaxy Gyr~' found by Boselli et al. (2008) for Virgo and 
the latter consistent with the estimate for Coma). 

An alternative method for estimating the infall rate can be 
achieved using the number of LTGs that are currently in the HI- 
poor phase divided by the time scale for complete hydrogen ex- 
haust: 0.5 Gyr for low-mass galaxies (Boselli et al. 2008) and 
somewhat longer for massive systems. This time scale is esti- 
mated assuming that a complete HI and H2 ablation will occur 
during the second passage through the cluster center. This is con- 
siderably shorter than ~ 2 + 1 Gyr (Bigiel et al. 2008): the time 
scale for complete consumption of the H2 due to star formation 
in intermediate HI deficiency galaxies, i.e, without an extended 
reservoir of atomic hydrogen that can replenish the inner HI fuel 
consumed by the star formation (Fumagalli et al. 2009). The de- 
rived infall rates adopting these figures are consistent with the 
previous estimates. 

Finally, assuming that all galaxies (LTGs+ETGs) were ac- 
creted by the Virgo cluster at a rate that remained approximately 
constant in time, and that massive galaxies require an additional 
1-2 Gyr after a stripping event to migrate from the blue cloud 
to the transition region and to the red sequence (Boselli et al. 
2008; Cortese & Hughes 2009; Gavazzi et al. 2010), we infer 
that the age of the Virgo cluster is roughly 2 Gyr, implying a 
cluster formation at z ~ 0.2, much later than z ~ 1 derived using 
similar arguments for the Coma cluster in Paper EH. Of course, 
these are just order-of-magnitude estimates of the accretion rate 
on the Virgo cluster. For example, cluster B (M49) will soon 
merge with cluster A (M87) to form an even more ralaxed and 
rich Virgo cluster, contributing with many galaxies that have 
been pre-processed in the group itself. The significant age dif- 
ference found between Virgo and Coma is consistent with the 
high contrast in the observed X-ray properties of the two clus- 
ters (Forman & Jones 1982), suggesting a marked difference in 
their evolutionary stages. 

Table 4. Build-up of the Virgo cluster from various estimators. 
The first number refers to objects with stellar mass > 10 7 M Q , 
the second to those with > 10 9 M (comparable with Coma) 





All 




Hl-poor 


Hl-rich 


N (> 10' M Q ) 


738 




160 


93 


t (Gyr) 


2 ± 1 




0.5 ± 0.4 


0.2 ± 0.1 


Infall rate (Gyr 1 ) 


369 ± 


184 


320 ± 190 


465 ± 230 


N (> 10" M s ) 


210 




73 


36 


t (Gyr) 


2 ± 1 




0.7 ± 0.4 


0.3 ±0.1 


Infall rate (Gyr 1 ) 


105 ± 


53 


104 ± 61 


120 ± 45 



8. Summary and conclusions 

The wealth of data provided by ALFALFA in the stretch of the 
Local Supercluster visible in spring, which contains the Virgo 
cluster, which were followed-up with Ha imaging observations, 
enabled us to analyze the relations between the HI gas, the stel- 
lar content and the rate at which the gas is currently fueling the 
formation of new stars in spiral galaxies. 

Given the sensitivity of ALFALFA and of Ha3 at the distance of 
Virgo, our analysis comprises galaxies with a stellar mass in the 
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Fig. 10. (Left panel) Relation between the specific star formation rate and the gas fraction in three classes of increasing Defni, 
excluding objects fainter than M» = 10 9 M Q . (Right panel) Same as in the left panel, but including fainter objects. The histograms 
drawn on top and to the right show the counts of galaxies with given HI gas fraction and SSFR, with the mean values given in three 
bins of Defni- 



range 10" <M, < 10 1L5 M , a gas mass 10" < M H , < 10 10 
M and a star formation rate in excess of 10~ 3 M yr~'. The main 
results of the present investigation can be summarized as fol- 
lows: 

- In addition to morphology segregation, the Virgo cluster rep- 
resents a significant perturbation on the Local Supercluster 
for the gas content and the star formation properties of its 
member galaxies of late-type, such that on average the HI 
content is reduced by a factor of four and the specific star 
formation rate by a factor of two. 

- LTGs in the Local Supercluster have progressively redder 
color with decreasing HI gas content, from normal defi- 
cient galaxies to objects that are undetected by ALFALFA. 
The gas-ablation mechanism also governs the star formation 
properties. 

- Highly HI deficient and undetected LTGs have ongoing nu- 
clear and circumnuclear (truncated) star formation, while HI- 
rich LTGs show star formation extended on the disk scale. 
We propose that the quenching mechanism responsible for 
these features is ram pressure. 

- The fundamental scaling relation (studied in Section [8]l be- 
tween the HI mass and the stellar mass in normal spiral 
galaxies is such that Mhi ~ M^ 5 . This implies that increas- 
ingly more massive galaxies have progressively lower HI gas 
fraction, consistent with the idea that they have transformed 
most of their gas into stars at higher z, in agreement with the 
downsizing scenario. 

- The scaling relation between Mhi and M„ followed by galax- 
ies in the Virgo cluster is significantly offset toward lower 
gas content, compared to that of Hi-normal galaxies. 

- Once the scaling law with stellar mass is removed by an- 
alyzing the correlation between the specific star formation 
and the gas fraction, it is found that a residual second-order 
correlation exists between these quantities. The effect of the 
environment is ten times stronger at producing the gas deple- 
tion than at quenching the star formation. 



- The present infall rate of ~ 400 galaxies per Gyr (of > 
10 7 M objects), inferred by the number of healthy spirals 
present in the Virgo cluster, is consistent with the number of 
galaxies of all types (of similar mass) if the infall process 
has been acting for approximately 2 Gyr, about three times 
shorter than on the Coma cluster (Paper III). This is consis- 
tent with the idea that Virgo is a young cluster caught in an 
early evolutionary stage. 
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